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Quantum phase transitions are often embodied by the critical behavior of purely quantum quan- 
tities such as entanglement or quantum fluctuations. In critical regions, we underline a general 
scaling relation between the entanglement entropy and one of the most fundamental and simplest 
measure of the quantum fluctuations, the Heisenberg uncertainty principle. Then, we show that 
the latter represents a sensitive probe of superradiant quantum phase transitions in standard mod- 
els of photons such as the Dicke Hamiltonian, which embodies an ensemble of two-level systems 
interacting with one quadrature of a single and uniform bosonic field. We derive exact results in 
the thermodynamic limit and for a finite number A'' of two-level systems: as a reminiscence of the 
entanglement properties between light and the two-level systems, the product Aa;Aj3 diverges at the 
quantum critical point as N^^^ . We generalize our results to the double quadrature Dicke model 
where the two quadratures of the bosonic field are now coupled to two independent sets of two- 
level systems. Our findings, which show that the entanglement properties between light and matter 
can be accessed through the Heisenberg uncertainty principle, can be tested using Bose-Einstein 
condensates in optical cavities and circuit quantum electrodynamics. 

PACS numbers: 



I. INTRODUCTION 

Quantum phase transitions (QPTs) occur at zero 
temperature and are triggered by quantum fluctuations 
[T]. They manifest themselves in a sudden change of 
quantum properties of the collective ground states. The 
appearance of a flnite order parameter or collective 
gapless excitation is often used as a signature of such 
a critical phenomenon. Recently, tools of quantum 
informations were also used in this perspective. For 
instance, the Von Neumann entanglement entropy 
[2] exhibits a critical behavior in many collective 
models undergoing a QPT [3]-[5]. Other examples of 
such quantities can appropriately locate the Quantum 
critical points (QCP) of some given systems. Among 
them, one can quote the fidelity (or overlap) [HI which 
directly quantifies the suddenness of the change of 
the ground state wave-function around the QCP, the 
global geometric entanglement [7] which evaluates its 
distance, in the Hilbert space, to the closest separable 
state, the logarithmic negativity ISJ, which measures the 
entanglement between non-complementary parts of the 
system, and more recently, the bipartite fluctuations 
which have been shown [S] to be directly related to the 
entanglement properties between subsystems of some 
one- or higher- dimensional fermionic models. 
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This article is aimed to demonstrate that the critical 
behavior of the Heisenberg principle (HP) can be used 
as a characterization of QPTs as well as entanglement 
properties, in particular in systems of photons. 

Since QPTs are driven by quantum fluctuations of the 
flelds that define the physical structure of the system, 
it is rather natural to focus on the simplest measure of 
those fluctuations, the product Aa;Ap. Secondly, and 
importantly, fluctuations appear in the calculation of the 
entanglement of many important physical models |10j . 

In the case of light-matter systems, one can derive 
the Von Neumann entropy of the reduced density ma- 
trix pr of the photon mode S = — Tr{pr log2(pr)}- The 
latter is obtained at zero temperature by tracing out 
the matter degrees of freedom in the density matrix 
p — \G){G\ (where \G) is the collective ground state), so 
that pr — Trj^g^^^gj.{|G) (G|}. Then, providing that the 
model can be described by a quadratic Hamiltonian of 
interacting bosonic fields llj, the entanglement entropy 
S will be directly related to the HP through: 

S = -Tr{p,.log2(pr)} 

= (^AxAp + log2 I^AxAp + i| 

- (^AxAp - log2 I^AxAp - i| , (1) 

with AxAp = (l/2)V(G|(a + at)2|G)(G| - (a-at)2|G) 
{h = 1). Since we consider that this equality constitutes 
an important relation between entanglement properties 
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of the system and fluctuations, for completeness, we shall 
derive it below, in Sec. [ll] From this equation, one even 
obtains the very simple relation: 



S log2{AxAp} if AxAp > 1. 



(2) 



Thus, a divergence of such a quantity would directly 
imply a logarithmic divergence of the entanglement 
entropy. The a-Renyi entropies, which can individuate 
a QCP [12], are also very sensitive to the divergence of 
the HP. Finally, the individuation of a QPT through 
the HP seems to be an important concept to push 
forward owing to the fact that such quantity is directly 
accessible experimentally, especially in the context of 
quantum optics. Incidentally, it is remarkable to notice 
all the efforts made so far to develop methods for 
determining the entanglement in many-body systems 
from the measure of physical observables (T3j . 

In the following, we will investigate two simple yet 
relevant models exhibiting QPTs. First, the HP will be 
investigated in the case of the celebrated Dicke model 
|14j . where a bosonic field interacts through one of 
its quadratures with a chain of N two-level systems. 
Then, the double quadrature Dicke model [TH], where 
the two quadratures of the bosonic field are coupled 
to two independent chains of atoms, will be analyzed 
thoroughly. This recent model provides a very natural 
framework for the investigation of the critical behavior 
of the product AxAp at the QCPs of light models, 
because the fiuctuations of the quadratures x and p, 
which both interact with two different sets of atoms, 
play an equivalent role and can both undergo a critical 
enhancement at the QPT. Finally, possible applications 
of those models using Bosc-Einstein condensates in 
optical cavities [TBHT5| or circuit QED P^H2^ wiU be 
briefly discussed. 



The paper is organized as follows. In Sec. [TTj we will 
review the method of calculating entanglement entropies 
in quadratic Hamiltonians of interacting bosonic fields in 
order to demonstrate Eqs. ([T]) and Then, in Sec. 
Ill we will explore the criticality of the photonic fluctu- 
ations in the Dicke model, both in the thermodynamical 
limit case and for a finite number N of two-level sys- 
tems. In Sec. |IV[ the fluctuations and the entanglement 
in the double quadrature Dicke model will be investi- 
gated, with a special emphasis on the point of double 
symmetry breaking where those quantities behave in a 
peculiar way. Finally, we conclude in Sec. |V] 



II. RELATION BETWEEN THE 
ENTANGLEMENT ENTROPY AND THE 
FLUCTUATIONS IN QUADRATIC MODELS 

We review here the derivation of the entanglement en- 
tropy and the a-Renyi entropies for the ground state of 
a quadratic Hamiltonian H(a, 62, 63, &«) expressed in 
terms of n different bosonic fields a, &2, ^3, ^n, by using 
the method introduced in the articles [11] . Note that the 
following method can also be used for fermionic modes 
[11| . First of all, it is well-known that the ground state 
of any quadratic bosonic Hamiltonian can be written as: 



\G) =AAe(*')^^*'|0)® |0)..|0), 



(3) 



where (3"^)^ = {a' ,bl,bl, ...,bl^), ^ is a symmetric (real 
or complex) nxn matrix , A/" is a normalization constant 
(for instance real), and |0) (g) |0)..|0) is the vacuum of 
the fields a, 62,^3, ■•■,&« : a\0) (g) |0)..|0) = and bi\0) 
|0)..|0) — for i — 2...n. We now derive the reduced 
density matrix of the photon mode a by tracing out the 
matter degrees of freedom &2, ^3, in the pure ground 
state density matrix p = \G){G\: 



Pr = TrjjiatterM = Trb,,..fc„{|G)(G|} 
r2 f f f f du2dv3...dvn 



(4) 



dv, 



u2\®{yz\...®{yn\e^'''^"^'^'\0)®\Q) 



\V2) ® \V3)--- ® Wn) 



r 



rn-2 



where for i — 2. 
for the i*'^ field: 



the states 



are coherent states [53] 
|0) (satisfying bi\iyi) = 



i^iji^i)). The index iq is such that Aug ^ (if such an in- 
dex does not exist, the problem becomes trivial). Then, 
we introduce the new variables v'^ = Vi (« 7^ «o) £^nd ~ 
2 Sr=2 -^li'^^- N^O'v^ since the integral of a gaussian is a 



gaussian, and using that J ^e'''*"' i-^' Vio)a('^jo U = 
\xfa+i J ^f^gy< \v^,)a{v^^a = for any num- 
ber X and 1", one can realize that pr is proportionnal 
to the product of several exponentials of a^, (a^)^ and 
a) a. One concludes !24J that p,. is the exponential of a 
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quadratic Hermitian form in the field a : 

Pr = cxp{— Kq — Kio^a — n2C? — ^2(0^)^}, (5) 

where kq and ki are real, and K2 a 'priori complex. One 
diagonalizes — {kq + Kia^a + K2a^ + K2(a^)^} by a Bogoli- 
ubov transformation to write : 



(6) 



where P = ua+ua^, with u and v such that — |wp — 1, 
(and u chosen to be real), where A is the pseudo-energy 
(not to be confused with the fluctuations Ax or Ap), 
and where is a constant. Then, since Tra(pr^) — 
Tra,fc2^..b„ (v4p) — {G\A\G) for any photonic operator A, 
it is easy to determine , A, u, and v. Inverting the 
Bogoliubov transformation leads to a = u*P — vP'^ and 
allows us to write: 



1 + e- 
1 - e- 



= ((G|a^a|G) + {G\aa}\G)f 

-4(G|(at)2|G)(G|a2|G) 
= 4{(AxAp)2 + C}, 



(7) 



where we have notably used that ■^_^-t^ = 1 since 
Tra(pr) = {G\G) = 1 and where we have introduced 
C = ((G|a2|G) - (G|(at)2|G))V4. Finally, the pseudo- 
energy reads : 



A = loe 



'2^(Aa;Ap)2 + C + 1' 
2v/(Aa;Ap)2 + C- 1 



(8) 



The derivation of the entanglement entropy becomes 
now straightforward; 

S = -Tra{Pr log2(Pr)} = -Trp{pr loga (pr)} 



k=0 



fro log (2) 

{^{AxApf + C + ^)log2 



(9) 



(AxAp^+C + \} 



{ViAxApy+C ~ -) log2 \ ^{AxApY + C 



where the states \k) are Fock states for the P operator. 

A priori^ the parameter C, is not always equal to 
[IS]. But on the other hand, it is always possible 
to come down to the case ^ = . In fact, if one 
starts from a quadratic Hamiltonian 'H(a, 62, bn), with 
(G|a2|G) = e^"l'\{G\a?\G)l (and 2</> ^ 0[7r]), by in- 
troducing a — ae"**^, the formally new Hamiltonian 
H(a, 62, 6n) = 'H{a,b2, ■■■,bn), has a ground state 
|G) = AAe(*')^-^*'|0) |0)..|0) = |G), where (($)t)^ = 
((o)t,6^,...,6j,), and where A = i:»(e-*^)yli:»(e-*^) with 



D(e *'^) the n x n diagonal matrix defined by £)(1, 1) = 
e-"t' and i:'(j,j) = 1 V2 < j < n. Then {G\a^\G) = 
|(G|a2|G)| and4(Ai)2(Ap)2 = -{G\[d+{a)^]''\G){G\[&- 
(a)t]2|G) = 4(A2;)2(Ap)2 + {{G\a^\G} - (G|(at)2|G))2. 
Finally, even if the introduction of a new photonic oper- 
ator a is needed, one can always convey to Eq. ([!]). Note 
that for the Hamiltonians studied in the present article, 
i.e. the standard Dicke Hamiltonian (See Eq. ( [l2| ) below) 
and the double quadrature Hamiltonian (See Eq. ( 18 ) be- 
low)), one directly has C = 0, with no need to introduce d. 
In the case where AccAp diverges, the Taylor expansion 
log2{Aa;Ap± 1/2} = log2{Aa;Ap} ± [2 log(2)Aa;Ap]-i + 
o[{AxAp)'^] implies Eq. ([2]) : S" - log2{Aa;Ap}, which 
provides a very simple relation between entanglement 
and fluctuations around the QCP. Equivalently, the a- 
Rcnyi entropies, defined as: 



Saipr) 



1 



1 



log2{Tr(p«)} 



(10) 



might also be written as: 



Sa{pr) 



a - log2{(l + 2Aa;Ap)" - (2Aa;Ap - 1)"} 



1 - a 



(11) 



Note that the entanglement entropy corresponds to 
Sa^x. In the case where AxAp diverges, the a- 
Renyi entropies also diverge as a logarithm : Sa{pr) ^ 
log^iAxAp). 

Below, to illustrate this relation between entanglement 
properties and the Heisenberg principle, we focus on stan- 
dard models of photons and compute directly the HP. 



III. THE DICKE MODEL 

We first focus on the Dicke Hamiltonian (DH), which 
describes the coupling of a single and uniform bosonic 
mode a of energy uj with N two-level systems with atomic 
splitting Wo (again, the Planck constant h has been fixed 
to unity for simplicity) [Hj : 

H = Lja}a + + ^={a + a)){J+ + J^), (12) 



where A is the atom-field coupling strength. The total an- 



gular momentum operators and J± read = '^z 
and J± = T^^=i^^±-i where cr^ and CTj_ are the usual 
Pauli matrices for the l*-^ pseudo-spin, so that the an- 
gular commutation relations are [J2, J±] = ±2J± and 
[J+, J_] = 2Jz. In the thermodynamic limit (A^ — )■ 00), 
this Hamiltonian undergoes a superradiant QPT for A — 
A"^^ = ^wwg /2 j2B]. When A < A'^'', the system is in a 
Normal Phase (NP), with a squeezed and non-degenerate 
vacuum, while the Superradiant Phase (SP) occurring 
for A > A'^'' is embodied by the appearance of a double 
degeneracy with atomic and photonic macroscopic co- 
herences. This QPT is associated with the breaking of 
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FIG. 1; (Color online) In red, Aa;Ap (HP) for the ground 
state(s) of the Dicke model (cf Eq. (12l) in the thermody- 
namic limit (A'^ — >■ oo), with the Planck constant h fixed to 
unity. As a reminiscence of the entanglement entropy between 
atoms and light [3], it diverges at A = A'^'' as |A — A'^'']"^''''. 
In dashed line, results given for the finite-size ground states 
of some exact diagonalizations. Since those are Schrodinger's 
cat like (restoring the broken symmetry), their HP diverge 
when A ^ A"^"^, in contrast to the thermodynamic case. Inset: 
AxAp at A = A'^'^ versus A'^ (the number of two-level systems), 
for = 10 to iV = 10^ it scales like N^'^ |3Sl|37j. Note that 
the limit of A'^ = 2.10^ ^'^Rb atoms coupled to an optical cav- 
ity has been achieved in BEC realizations of the Dicke model 

Ezl- 



the parity operator n which accounts for the parity of 
the total number of excitation quanta: 11 — exjp{iTT Nf^^c) 
with Nexc — cL^a + Jz + N/2. Moreover, through the 
Holstein-PrimakofF transformation [27], which allows to 
write the angular operator in terms of a bosonic field 
b: J+ = 6t(iv _ &t;,)i/2 ^ j_ ^ (TV - b^by/'^b and 
Jz = b^b — N/2, the DH can be proved to be equiv- 
alent to Hjy^oo — A+e^_e+ -I- A_e!Le_ -|- Eg- 

Eg is the fundamental energy, and the normal eigen- 
frequencies (gaps) are such that 2Aj. — (wq//^)^ + cj^ ± 



v/((Wm? 



w2)2 + WX^LuuJon, with /.t = 1 if A < A^'' 
The bosonic operators (the po- 
laritons) e± are linear combinations of the original fields: 



and M = if A > A' 



-j§tf[{A^-wo){b^ + (A_ +c^o) (&-/?)] I (13) 

_ 1 f sin(7) 
21 



[(A+ - w)(at +a) + (A+ + Lu){a + a)] 



+ 



COS (7) 



[(A+ - Qo){b^ -f3) + (A+ + Cjo){b - 13)]}, {U) 



where ujq — + l)/(2^), a = eyjN{\ — ^'^)X/uj is 

the photonic coherence, (3 — t\jN(l — /i)/2 the elec- 
tronic coherence, with e = in the NP and e = ±1 
in the SP. The mixing angle 7 is such that tan(27) = 
(4A ^loloq^^/'^) / (wq — ^i^ijJ^). The coherences are zero in 



the NP while they can be either positive or negative in the 
SP, resulting in the double degeneracy of the eigenspec- 
trum. In both cases, the ground state(s) \G)^, defined 
for a given set of coherences are the ones of a double har- 
monic oscillator (shifted in the SP), and satisfy e_|_|G)e = 
e„|G)e — [28 . Besides, since 11 — exp(i7r{a^a + b^b}), 
then nall^ — —a and 11611^ = —b. Consequently, in the 
SP, n|G)± = |G)ip , i.e. the symmetry 11 is broken : in 
contrast to the NP, the ground states \G)± are no longer 
eigenstates of the parity operator 11. After inverting the 
last polaritonic relations, one gets, for each ground state 



AxAp - 



V 



2A_ 



2A+ 



}{ 



cos(7)^ A_ 



sin(7)2 A_|_ 
2 



},(15) 



which is plotted in Fig. [l](in red). For A — > A^*", Ai ^ 
I A - A'='~|i/2 [26] and the HP diverges like |A - A'^^^^^'^- 
Thus, in the thermodynamic limit, where the DH in Eq. 
(12) is quadratic (see Hn^oo above), the criticality of 
the ground state entanglement [3] at the QCP is indeed 
involved by the divergence of the HP, as proved by Eq. 
([2]). It is also important to compare the criticality given 
by the mean-field with some finite size results. First, 
in the limit A/w — )■ (X) (the so-called ultrastrong cou- 
pling limit [29U31) ). an TV*'' order perturbative theory 
allows us to prove that the two first eigenstates |5'g) and 
I^Pb), have their energies separated by an exponentially 
small splitting, and are linear superpositions of the states 
la^) 1^/2)2; and l — ap) \ — N/2)x- Here, \±aF) are coher- 
ent states for the photonic part: a\ ± ap) = ±a;i?| ± ap) 
with ittF = ±y/~N\/uj [31 [33]. The states | ± 7V/2)^ 
are the two maximally polarized Dicke states in the x- 
direction (the direction of the coupling) : they satisfy 
U±N/2) = (1/2)(J+ + J_)|±Ar/2) =±iV/2|±Ar/2) = 
±A^/2H^j^|±)j where each local pseudo-spin state |±)j 



satisfies al\±)i = ±|±' 



3- 



The light-matter coupling 



is so important that each pseudo-spin is polarized in 
the direction of the coupling ( the x-direction in this 
paper). Moreover, since H|ai?) 1^^/2)2: = e"" '^I'^-f) 



lowing cat's states wave-functions [3] 
that restore the broken symmetry H 



- N/2).x, the fol- 
are the only ones 



^75{i«F>iA^/2). + (-in- 

^j-^{\ap)\N/2)x-{-ir\ 



ap) I - 



N/2)x] 
iV/2)4l6) 



Then, from the expression of \^g) in Eq- ( 16 ) one iden 



tifies AxAp ~ {a% + 1/4) ^ y/NX/u which asymp- 
totically matches the finite size curves of Fig. [T] for 
\/lj 3> 1. In fact, Y^g) is a symmetric superposition 
of the two thermodynamical vacua \G)± which are the 
ground states of a double harmonic oscillator shifted 
around some macroscopic coherences. Since those coher- 
ences get infinitely far from each other when increasing 
both the number of atoms N and the coupling A, the fluc- 
tuations of such a superposition diverge also in this limit, 
contrary to the thermodynamic limit result. Finally, we 



5 



must evaluate the fluctuations of the finite-size ground 
states at A = A'^'' to see whether it diverges, or not. The 
scaling hypothesis which relates the exponent of N and 
the power of |A — A'^''| in the finite-size developments of 
every physical observables at the QCP allows us to 
prove that : 



AxAp-iV^/s for A^^A'^ 



(17) 



in quantitative agreement with the numerical simula- 
tions, confirming the divergence of the HP when TV — >■ cx). 

Actually, in the standard DH, the critical scaling of 
the HP comes from the criticality of Ax because the 
quadrature x is the one that interacts with the two- 
level systems. But what happens if the two quadratures 

- a^) are 



X = (l/(V2w))(a -I- aT) and p = 
coupled to two different chains of atoms? 
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IV. THE DOUBLE QUADRATURE DICKE 
MODEL 



Let us consider the double quadrature Dicke Hamilto- 
nian which has recently been introduced ^ISj: 




(18) 



i^=(a-a^)J^, 



where the chain of two level systems labeled by C (resp. 
/), has an atomic transition frequency (resp. ujj ) 
and is coupled to the quadrature a + (resp. i[a — a^) 
via the coupling strength constant Ac (resp. A/). 

Here, two independent symmetries transformations 7/ 
and Tc are conserved and are defined via the following 
operations: 

Ti 



(a + a' , i{a — a' ), , J^) {a + a' , —i{a — a' ), , —J^) 

{a + a\i{a - a^), jf , 4) ^ {-a - i(a - a^), - , 4) > and A/ < Af , the lower energy gap Ai vanishes as 



FIG. 2: (Color online) Bottom left panel: two-dimensional 
phase diagra m o f the double quadrature Dicke model [see 
Hamiltonian ( 18 1] . Other panels: fluctuations Ax (black 
dashed dotted line), Ap (blue dashed line) and Heisenberg 
Principle AxAp (solid red) in the resonant case (cJcau = i-^c ~ 
lJ^ = 1), in the thermodynamic limit {Nc,Ni — >■ oo), and 
with the Planck constant h fixed to unity. The resul ts are 
plotted with respect to the radial coupling ^ A^ -I- A|, with 
\i/\c = est — tan(6'), for several polar angles 9. At the 
QCPs, the product AxAp shows either a divergence scaling 

as Ai, (where Ai is the lower gap) which is reminiscent 
of the one in the standard Dicke model, or a local maximum 
(at the point of double symmetry breaking, for 6 = 7r/4). 

The gaps Ai < A2 < A3, the polaritons ei, 62 and 63, 
and the fundamental energy Eq are obtained by diago- 
nalizing the associated Bogoliubov matrix |15) . 

We first study the case Ac > A/ , which corresponds to 
9 < 7r/4 where tan(6') = A//Ac (see figure J2]). For Ac ^• 



where a^a, and remain unchanged. 7/ can be 
viewed as the time reversal symmetry |15) . Again, 
by an Holstein-Primakoff transformation, one introduces 
the bosonic fields be and bj defined via the relations: 



I Ac 



Apl^/^, and the lower polariton ei reads: 



ei 



Ho 



^?Ai 



Jl = [Nk - blbk 



and 

One can then show 



4 = bliN, - blb.y/^ 

= blbk - Nk/2 {k e {I,C}). 
that 7c and 7/ gets broken w hen Ac and Aj are i n- 
creased above Ac — -y/wcau wJi,/2 and XY — \J ^cav^'il'^ 
[15j . This gives rise to four different quantum phases in 
the thermodynamic limit {Nc-,Ni — > 00), separated by 
two orthogonal critical lines of equation Ac = Ac and 
A/ = Xf^ as shown in the bottom left panel of Fig. [2] 
One has either one (normal phase), two (real/imaginary 
superradiant phase) or four (double superradiant phase) 
degenerate and coherent vacua whose eigenfunctions are 
the ones of a triple harmonic oscillator : 



2i- 



^'l^cav ^ 4A| 

o^A/Ai 



4Xj 



bi-il- 



x-(l- 
c^?Ai 



Ai 



■ be 



4Xj 



(1 



r- 



cav ^c 



2i 



y/^A/Ai 

^j^cav ^ 4Aj 



^1}, (20) 



where TVei = 2{Ai(- 



-4\j 



'-)} 



1/2 



Eq. (13l). In particular, the HP diverges for Ac ^ ^c 



H 



N- 



Aielei + A2ele2 + Aae^ea + . (19) ^cr 



Those polaritonic coefficients exhibit the same diver- 
gence as the lower polariton e_ of the standard DH (see 

^1^. Interest- 



(and stiU Xi < Xf) as 1/V Ai = lAc-AJ?]- 
ingly, we realize by examining the expression ( 20 ) , that 

the polaritonic divergence disappears when Ac — 
\i<\Y 



XY simultaneously, which occurs only 
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FIG. 3: (Color online) Entanglement entropy between the 
photonic field and the two matter fields S = — Tr(pr logj Pr) 
(where pr = Tibc,bi{\G){G\}) , with respect to ^^/X^ + Xj, 
at resonance (ujcav = i^c = "^i) along the radial lines 
A//Ac = tan(e') for 9 = 57r/16 (left), and 6^ = 7r/4 (right). 
Results shown for TV = Nc ~ Nj, in the thermodynamic limit 
and in the finite size case. For N = 00 and 9 ^ n /A, S diverges 
at the QCP hke (-l/4)log2 |Afe - X'tf] {k = C,I), i.e analo- 
gously to the standard Dicke model [1]. For 9 = tt/A such a 
divergence disappears thanks to a compensation between the 
two simultaneous QPT (see Eq. (|20[ ) and ( |21[ )). Insets: same 
quantity at the QCP, versus iV (up to N=512). For 61 / 7r/4, 
S ~ log2(iV)/6, as in the standard DH [36l; for e = 7r/4, S 

does not diverge anymore: S{N) ^^Jl^ 5'<x,(Ac, A") — 1-29, 
providing a situation where a second order QPT admits a 
finite entanglement at its QCP. 

at the point of double symmetry breaking (at the cross- 
ing of the hnes Ac = A^ and A/ = ^Y)- ^^lis point, 
corresponding to 6* = 7r/4 if — w^, the two matter 
modes play a symmetric role and one has: 



ei 



ibi 



(21) 



There, the fluctuations Aa; and Ap do no longer di- 
verge and AxAp = 1/2 + (1 + 4(w[!,/wca,;)^)"^/^- As 
shown in Fig. [2] this value is a local maximum of the HP 
along the line Ac = A/. 

More generally, at every other QCPs of the two- 
dimensional phase diagram, where 7i or 7c are indi- 
vidually broken, the fluctuations of the quadrature in- 
volved in the transition diverge as 1/\/a7, while the 
other quadrature fluctuations remained bounded. Thus, 
at those single QCPs, one has either Ax ~ |Ac — Ag" |~^^* 
and Ap < 00 for Ac -> A^ or Ap \Xi - Af]"^/^ and 
Aa; < 00 for A/ — ■ Those scalings imply the critical 
behavior of the HP, which is the appropriate measure to 
detect all the QPTs in this model. 

Finally, we could illustrate our initial statement about 
the equivalence of the measures of the HP and the en- 
tanglement entropy S, by showing the behavior of this 
latter in this double quadrature model (See Fig. [s]). In 
order to compare the entanglement S in the thermody- 
namic limit to some finite-size situation, for which the 



ground states are cat's states in the superradiant phases, 
one must add to the expression Eq. ([I]) , a term account- 
ing for the degeneracy g], and equal to I (resp. 2) in the 
real or imaginary (resp. double) superradiant phases. 
That is why S saturates at 2 when Ac, A/ 3> ujcav 
As expected, for 6 7^ 7r/4, at the single QCPs (Ac or 
XT), the divergence of the HP ( due to the divergence 

involves the 



(20)) 



of the polaritonic coefficients, see Eq 
following scaling of the entanglement entropy in the ther- 
modynamic limit: 



5^(-l/4)log2|Afc-A^ 



(22) 



with k € {C, /}. This is reminiscent of the standard 
Dicke model [4]. Moreover, in the finite-size situation, 
the entanglement entropy scales as 5 ~ (1/6) log2(Af) 
(where N = Nq = Nj) at the single QCPs, as an other 
reminiscence of the standard Dicke model [36]. In the 
left panel of fig. |3j we show the plot of the entanglement 
S along 9 = 57r/16, both in the thermodynamical and 
finite-size cases. We clearly observe the two consecutive 
critical enhancement of S at the two consecutive QCPs. 
On the other hand, at the point of double symme- 
try breaking (corresponding to — tt/A at resonance 
jj) , the polaritonic coefficients do not 



diverge anymore (see Eq. (21)). Consequently, the en- 
tanglement S stays bounded in the thermodynamic limit, 
and apart from the discontinuity equal to 2 immediately 
after the double critical point (Ac , XY) ( due to the ap- 
pearance of the four fold degeneracy), its value, given by 
Eq. ([1]), reads at resonance : 



{1 + 75} 



V5 



l0g2(l 



1.29. 



(23) 



Correspondingly, for < 00 at this point, the physi- 
cal quantities admit a standard 1 /N finite-size expansion, 
and SiN) ^ 5oo(Ag:,Af ) < 00. Thus, while it under- 
goes a second order QPT, the entanglement of the sys- 
tem does not diverge when the number of pseudo-spins 
tends to infinity, which is somehow unusual, but in per- 
fect agreement with the behavior of the HP. 



V. CONCLUSION 

To summarize, in this work we have exemplified 
the enhancement of fiuctuations at superradiant QPTs, 
through the HP. By exhibiting a general relation valid for 
any quadratic bosonic Hamiltonian (see Eq. ([I])), we have 
shown that the HP is indeed connected to the logarithmic 
enhancement of the entanglement entropy, while being 
certainly easier to measure, since it does not require the 
full tomographic determination of the density matrix |38j , 
but just the measurements of the variance of the two or- 
thogonal field quadratures. By the way, we would like to 
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point out that the two models presented above could be 
physically implemented either with atomic Bose-Einstein 
condensates in optical cavities [T71 [31] or in circuit QED 
[T51 [351 HO] . For the latter proposal, the coupling to the 
quadrature {a + a^) is provided by the capacitive cou- 
pling of the quantized charge of a Josephson atom and 
the quantum voltage of a resonator [3UI, while the cou- 



pling to the quadrature i{a — a^) is made thanks to the 
inductive coupling which connects the resonator current 
to the flux of the qubit [3D] . 
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with C. Ciuti and J. Vidal. This work was supported 
by NSF under the grant DMR-0803200 and also by the 
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